NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



TECHNICAL NOTE 2600 


STRESSES AND DEFORMATIONS IN WINGS 

SUBJECTED TO TORSION 

By B. F. Ruffner and Eloise Hout 

* 

Oregon State College 




AfiDC 

IWRV 

AFL gglf 


Washington 
February 1952 




' • ( 

I 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

TECHNICAL NOTE 2600 

STRESSES AND DEFORMATIONS IN WINGS 
SUBJECTED TO TORSION 
By B. F. Ruffner and Eloise Hout 

SUMMARY 


tech library kafb, nm 



DDLSS4S 


Basic equations of Karman and Chien (given in "Torsion with Vari- 
able Twist," Jour. Aero. Sci. , vol. 13, no. 10, Oct. 191+6, pp. 503-510 ) 
are solved by representing the shape of a torsion box by means of a 
Fourier series. The coefficients of the series are determined by con- 
ventional Methods. Angles of twist, longitudinal stresses, and shear 
stresses are determined in terms of the series coefficients. The 
method is applied to the calculation of angles of twist and stresses in 
torsion boxes of rectangular, elliptical, and airfoil cross section. 

Results obtained for angles of twist and normal stresses are in 
good agreement with results of Karmen and Chien except at sharp comers. 
Results obtained for shear stresses indicate the necessity for use of 
a large number of terms of the series for satisfactory accuracy. 


INTRODUCTION 


Experimental and theoretical investigations related to torsional 
stresses and deformations are given in references 1 to 18. In refer- 
ence 1, Karman and Chien have developed equations applicable to the 
problem of restrained torsion of tubes of arbitrary constant cross 
section loaded with a couple. Their solution of these equations was 
given for rectangles. In general, however, the method of solution of. 
the equations given by them is not easily accomplished. It is the 
purpose of the present report to give the' theoiy of an approximate 
method of solution of these equations and the results of the application 
of this method to tubes of rectangular, elliptical, and airfoil cross 
sections. This work was conducted at the Oregon State College under 
the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. " 
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SYMBOLS 


In figures 1 to 5 are shown coordinate systems adopted for this 
report . 

x coordinate along axis of tube 

y» or dinat e of airfoil measured from chord 

s distance along periphery of cross section from some 

arbitrary point 0 1 

r no rmal distance from shear center of section to tangent to 

wall at s 

r Q distance from leading edge to shear center 


t 



x 


u 

0 

E 

V- 

G 


I = 


thickness of wall of tube 

normal stress in skin in direction parallel to x-axis 
normal stress in skin in direction perpendicular to cr x 
shear stress in skin in direction of x- and s-axes 
normal strain in skin in direction of x-axis 

displacement of any point in skin in x-direction 
rate of twist of tube per unit length in x-direction 
Young's modulus 
Poisson's ratio 

I 

shear modulus 



ds 
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Ps 

t = / y' ds 

u 0 

k 2 . 

2 

a, b dimensions of rectangles and ellipses (see figs. 3 and k) 

5,- r) coordinates of point on ellipse (see- fig. I 4 .) 

0 = sin - -'- ~ (for ellipse) 
b 

n^ m, i positive integers 

V B n> C n> V B m> C m series coefficients 


v n 


coefficient 


enclosed area of cross section | i’fp r ds J 


total length of cross section X Q) ds 


M 


twisting moment 


THEORY FOR TORSION TUBE OF CONSTANT CROSS SECTION 

Assumptions 


The development of the equations of Kannan and Chien are based on 
the following assumptions: 

(1) Torsion tube is of constant cross section of arbitrary- shape. 

(2) Tube is loaded with constant torsional couple acting about an 
axis perpendicular to the cross section. 

(3) Bending. stiff ness of thin walls is negligible. 
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(I4.) Displacement of any point on wall of the' tube is composed of 
a displacement due to rigid rotation of cross section in plane of cross 
section plus a displacement due to warping. The latter displacement is 
in a direction parallel to the axis of the cylindrical tube. 

(5) The torsional deformations can be assumed to be independent 
of deformations due to bending and shear loads. This implies that the 
principle of superposition is applicable. 

(6) Because of the presence of bulkheads, the strain in the circum- 
ferential direction is negligibly small. 


General Theory 

Basic equations.- On the basis of the above assumptions, Karman 

and Chien have shown that the following equations are applicable 
(reference 1, equations 7j 9, 11, and 12):. 


2 3^u + 3^u + Q dr _ 

1 ~ t 1, 3x^ a s 2 ds 


( 1 ) 



( 2 ) 



ds + 10 



(2a) 


from which, for constant wall thickness. 


1 3 2 u + 3^u + t dr 
k 2 3x 2 3s 2 1 ds 


dr . ds + JL dr = 0 
ds GI ds 


(3) 


The general problem is to find a function u = u(x, s) that will 
satisfy equation (3) and the boundary conditions. In ary particular 
problem r = r(s) is known from the geometry of the cross section. 

This may, however, not easily be expressed analytically. 

An approximate solution.- For the approximate solution, let the 

function r = r(s) be given by a trigonometric series of the following 
form: 
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f-c„' + 




COS 


2nns 


■ (10 


This series is periodic in the interval s = l. The coefficients may- 
be determined by a number of methods to be discussed later. 

Differentiation of equation (1;) with respect to s gives 


— = -S~ 2nn C n ! sin - ri — s - 
ds — n l 


( 5 ) 


Let C n = -2nn C n '. Then equation (3) becomes 


dr V _ . 2nns 

li = Z_ C n sin — 


( 6 ) 


A particular solution of equation (3) may be shown (reference l) 

to be 


^(s) 


JO_ s 
2GAt V 





The solution of equation (3) may then be written in the form 

. u = ^(s) + u 2 (x,s) (8) 

If equations (8) and (7) are substituted in equation {3), the 
following is obtained: - 


_1 5^2 + 



+ t dr - 

I ds 


u c 


dr 

ds 



( 9 ) 


If the condition that the end x = 0 is restrained, then at 
x = 0, u = 0 and u^ = -u 2 . 
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If at the other end, say x = <*> } the end is unrestrained, then 


ct x = 0, or since 


E 


E Ou 


J x p 

1 - ti 2 


_ u. 2 ax 


1 - p 


then at x = » 


3u 


3u2 


= 0 


3x 3x 

A function U 2 satisfying these boundary conditions may be written 


as 


u 2 


- -H 


-\ n kx . 2nns 
sin — ? — 


(10) 


The coefficients A n may be determined by the condition that at x = 0,' 
U 1 3 -u 2^ or 


sin 


2nns _ If l 


2GAt 


>S 

,2 I 1 


s _ l_ | /r\ (j/s.N 

l 2aJ 0 \ j / W 


Let 


(H) 


= 


2GAt 


“n . '( A n) 

Equation (11 ) may then be written 


( 12 ) 


% 


sin 


2nns 


f -4J 0 (?) if) <«) 


Equation (13) may be used to determine the coefficients ^ and conse- 
quently A^ when the geometiy of the cross section is known. 
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x 


The solution of equation (3) satisfying conditions at x = 0 and 
00 may then be written as 


u = u-^s) + UgCxjs). = ^ A n sin ^ A n e' 


,-X n fac 


sm 


2mts 


or 


-14- ^ nkX ) sin S25 (lit) 


Since equation (ll;) satisfies the condition at . x = 0 and x = °°, it 
then is necessary to determine whether this is a solution that is appli- 
cable elsewhere . 

The particular solution u-j_ satisfies equation (3) so that it 
remains to determine whether U2 satisfies equation (9). Now from 
equation (10) 


_1_ 

k 2 Ox 2 


- A n 2 A n e ^ nkX sin 


2nns 


and 



z 


lpi 2 n 2 


V 


-X n kx 


sin 


2nns 

l 


and equation (9) may, by use of the above relations and equation (10), 
be written 


S in + 21 lUJ!_ A n e - ^n 

z 

■^(Z A n e " Xn ' 


“^n^ sin 2nns 


- I 


~ Xnkx sin issa 


t dr 
I ds 


^ . ds = 0 
ds 


( 15 ) 
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Since dr/ds = y G n sin 2nns/l, the integral in the above equa- 

tion may be evaluated. Since 

/(l : o a W x “ ta 

then equation ( 15 >) may be written 


^ Xn ^n e 


-l n kx 


sin 



2 n 2 . -X n kx . 2nns 
A n e n sin — - — 



sm 


2nns) 


Yi w 


-X m kx' 


= 0 


( 16 ) 


Equation (16) would be satisfied if every term of the summation is 
zero. Let any one term be denoted by the subscript ' i. Then 


2 -7u-kx . 2nis lpt 2 i 2 . -Xikx . 2nis 
A ap . 1 si n + - — - — A_-e sm — - 


- Xi A^e 


sm 


Af 6 


| °i 


This may be solved for X^ as follows: 






~(^m Xj_)kx 


(17) 


In the development of the solution the ' s were treated as 
constants. However, equation (17) shows that if the X^'s are constant, 
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the basic equations are not satisfied exactly. The question then arises 
as to whether a value of X^ for use in equation (10) can give a solu- 
tion that approximates an exact solution to the desired degree of 
accuracy. Since any constant value of each X^ will satisfy the con- 
ditions at x = 0 and at x = <*, it would appear that a constant 

can be found that will give reasonably accurate solutions in the 
region x > 0 to x ^ 

As x approaches zero, the value of X^ obtained from equation (17) 
approaches 


X ± 

x o 




( 18 ) 


Other values of X^ 
tinguish from the above. 


were tried and were, denoted by X^ to dis- 

-*-00 

These were taken as 



In the subsequent calculations to be discussed later,- both of these 
sets of values were computed and used to compute the angle of twist and 
stresses. Good agreement with an exact solution was obtained for rec- 
tangular sections by use of X^ obtained from equation (l8). 

If satisfactory values for X^ are obtained, then, since the A^'s 

are obtainable from equations (13) and (ll), the displacement function u 
(equation (lij.)) is obtainable. 

The stresses ct x , a g , and x and the angle of twist 0 may be 
given in terms of the displacement function as follows (reference l): 


CT x = 


_E 3u 

2 3x 


1 - \i‘ 


( 20 ) 


o s = pa x 


( 21 ) 
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e = JL + u 3E . ds 
IG I ,r ds 


( 22 ) 


t = G^ + roj 


(23) 


If equation (U4) is taken to represent the displacement, then the 
above equations become 


°x ^ p 3 ^n A n e 

1 - \ir — 


-X n kx 2nns 
sm — - — 


CT s = H W 

1 ~ V- £ ~ 


-Xnkx _ - 2nns 
sin — = — 




v — ( -X n kx\ • 2nns v / V " n - 2nns\ 

/ S “ ~ Ml'n” 1 — ' 

f]T nA n (l - e'^ 1 ™) cos 2ffiS + re 


(20a) 


(21a) 


ds (22a) 


(23a) 


Assuming equation (ll;) is a satisfactory approximation of the 
displacement function, then it is apparent that the stresses and angle 
of twist may be expressed, in terms of the coefficients A n , C n , and 
X n which are obtainable from the geometry of the given cross section. 

The problem is then resolved to that of determination of these coef- 
ficients. There are several methods available for determining these. 


Rectangular Cross Sections 

In figure 3 is shown the notation adopted for rectangular cross 
sections. The shape of the rectangle may be expressed nondimensionally 
by the parameter 


P = 


(210 


a + b 
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The coefficients C n ' for n / 0 (equation i|) may be obtained by 
evaluation of the integrals 


C ' 
n 


0 


r 2nns j 

— cos — - — d 

v L 


(!) 


( 25 ) 


The length of the periphery of a rectangle is 

l = U(a + b) 

The values of r are 

0 ^ s a 


a s "= a + 2b 
a + 2b < s ^ 3a + 2b 
3a + 2b ^ s ^ 3a' + k*> 


r = b 
r = a 
r = b 


r = a 
r =» b 


(26) 


3a + i;b ^ s "= Ija + lib 
The ratio v/l may be expressed in terms of p as follows: 


? - ! - i (1 - e) 


r a _ 8 
l “ l h 


£ = a „ £ 
III 1 


r 

i 


b 

l 


= io - ?) 


!*°*£ 

£ >£< i _ £ 

I - h~ 2 h 

£>i_£<± + £ 

1 ~ 2 k ~ 2 I4. 

S;>1 + £<1_£ 

12 h~ h 


£ > 1 - £ < 1 
l - ), = * 


00 




The integral of equation (2p) may then be evaluated to give 


, = 1 ~ 2P 


C , = ± ~ ‘-y sin 

n nn 2- 


nng 


(27) 


( 28 ) 


when n is even. When n is odd, C^' - 0. 
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The coefficients B n of equation (13) may be obtained by evalu- 
ation of the integrals 



These coefficients are then 


{ 2 } - 1) 3in nng 

" <«)*( P - f> 2 ) 2 


(30) 


By use of the relation C n = -2nrtC n 1 , 


C n = -2nnC n 1 = -2(1 - 2p) sin 


(31) 


The angle of twist may now be determined. For a tube of circular 
cross section of radius b and wall thickness t Q , the angle of twist 
is given by 


0 - JL 

circle j q q 


where 


*o 



ds = 2nb^t Q 


If the angle of twist of a rectangle is divided by the angle of 
twist of a circle whose diameter is 2b and whose skin cross-section 
area is equal to the skin cross-section area of a rectangle. 


6 = 9 

6 circle " 


JL .+ t 

IG I 



A. 


, n (l - e^) sin 


(r 


C n sin 2 ™^ ds 


1 ) 
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and from equation (12) 



Also, since 

I 0 = 2nb3t o 

l 2 = l6(a + b) 2 
A = iiab 

I = i;abt(a + b) 
2nbt 0 = UCa + b)t 


the above equation becomes 





The integral in the expression may be determined since terms only of the 
form sin 2 d^|j haye values other than zero. After integration 


the angle of twist may be written 


M/I 0 <5 


1 + 


HP £ B n o n (i - M 


]3y use of equation (2U) this may be written 


9 

H/l 0 G 




( 32 ) 



NACA TN 2600 


lit 


In a similar fashion it may he shown that the normal stress is 


°x = 


E MZk 


2 UAbtG 


( 2x n b ) B n e 


-(2X n b) 


kx 

2b S i n 2ima. 


(33) 


If it is assumed that A 0 = nb , then the above equation becomes 


°x = 


Ek M A-o* o 1 


-( 2 v>) 


kx 

2b . 2nns 


1 - m 2 *0*0° lO-tb 


x: (Mv ' n ' 2b 


A nondimensional stress cr x ' may be defined as 


CT„‘ = O-, 


(l - H 2 )A 0 t Q B 
EkM 


(3U) 


For the rectangle with l;(a + b)t = 2nbt Q , equation (3l|)' becomes 


i _ 


X 2p(l - p) 


1 V~ fn, h \r> “( 2X n b )^ . 2nns 
y (2k n b)B n e sxn — 


(35) 


A shear stress may be compared with the shear stress in a circular tube 
of radius b of the same skin cross-sectional area loaded with the same 
torsional moment. For the circular tube 


Tn = 


M 


2 ^-o b o 


(36) 


For the general case, from equation- (23), 
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Using the displacement function, equation (ll;), and equation (12), 
equation ( 23 ) may be written 



Then, using equation (36), the following relation may be written 


T 




nB. 


n 


1 - 



2nns 

l 


> + 


U /r\ 6 

1 - 


(37) 


Elliptic Cross Sections 

In order to evaluate the series coefficients for elliptic cross 
sections, it is necessary to express r as a function of s. In 
figure 1; is shown the coordinate system adopted for the ellipses. With 
this notation. 


where 



dt 


k' 



(38) 


(39) 


It may also be shown that 



r 


m 
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Therefore 


dr (k') 2 a 

Ik 

L J ) 2 

ds b 

1 - 0c') 2 |) 2 _ 


(Ul) 


Let 

sin 0 = 1 (1*2) 

D 

Then equations (38), (1*0 ) , and (1*1) may be written as 


\jl ~ (k 1 ) 2 sin 2 0 • d0 


^1 - (k 1 sin 2 0 


dr _ (k 1 ) 2 a sin 0 cos 0 
ds b r- 9 

|_1 - (k 1 ) 2 sin 2 0| 

The integral in equation (38a) will be recognized as the elliptic 
integral of the second' kind (reference 2) E(a,0). Tables of values of 

this integral are given in terms of the parameter 0 and a, where 

a = sin -1 k' (1*3) 

The total length of the periphery of the ellipse l is given by 


(38a) 

(l*Oa) 

‘ G*la) 




(k 1 ) 2 sin^ 0 • d0 


l = Ub 


(Uli) 
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Since the parameter a defines the ratio of minor to major axis 
of an ellipse (equations (U3) and (39)), this will be constant for any 
particular ellipse. This value is chosen first. Values of 0 are 
then assumed, from which £ and r) may be immediately computed. 
Values of s/b are then obtainable from the tabulated values of the 
elliptic functions. It is then possible to determine s/l, r/Z, and 
dr/ds corresponding to any assumed value of 0 or £ . 

The coefficients C n for the ellipses are then determined by the 
l equation 


n _ o dr 2nns 

C = c. — SXn 

n n ds Z 

* UU 


-(f) 


0*5) 


The integrals in equation (U5) may he evaluated by several methods. 

i 

In order to evaluate the coefficients B , equation (13) is used: 
For the ellipse this may be reduced to 


z 


L sin ^ = 2. - JL 
n l l 2n 


These coefficients are then given by 


0 * 6 ) 




(hi) 


For the elliptic thin-walled tube. 


I = t r c 


ds = l*ba 2 F^|,k') 


m 


where F^,k’j is an elliptic integral of the first kind (reference 2), 
Also, 

A = nab (li9) 

and 


l = lib 


E (f’ k ') 


(5o) 
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M/I 0 G 




— 


-OMg 

1 - e 

h B n C n' ” 


t 


( 58 ) 


« -iiy 

- 8A Z- 


(2X n b)B n e' (2inb) ^ .in 2Sf 


( 59 ) 



RESULTS OF CALCULATIONS 
Rectangular Boxes 


Calculations for angles of twist and stresses were made for rec- 
tangular boxes with cross sections defined by p = 0.05, 0.10, 0.15, 
0.30, and O.ijO. Computations of the coefficients C n and B n were 

first made by application of equations (30) and (3l). Tabulated values 
of these are given in tables 1 to 6. 

Values of A,- and \i were computed by use of equations (18) 

- L o - L eo 

and (19) for values of p = 0.05, 0.10, and 0.15* For p = 0.20, 0.30, 
and 0.1*0, only An values were confuted. These values are listed in 
x o 

the tables in the nondimensional form (2\ n b) 0 . The angles of twist 

for rectangles with p = 0.0p, 0.10, and 0.15 are listed in table 7- 
These were computed from equation (32).' Values were obtained based on 
values of both X Q and X m . For rectangles with p = 0.20, 0.30, and 

0.1*0, angles of twist were computed using only X Q . In table 7 are also 

given values of 9/(M/lG) .based upon a value of Poisson's ratio of 0.3, 
giving k = 0.592. This is the ratio of the unit angle of twist to 
that given by the theory which neglects the effect of warping of the 
cross section. 
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After this shear center is located, the airfoil is plotted and 
values of r/Z against s/Z may be measured, tabulated, and plotted. 
The coefficients C 0 ' and C n ' are then determined from the equations 


and 



Equation (13) is used to evaluate the coefficients 
coefficient may be written 


Any one 



Values of 
the form 


X ^ were computed from equation (18). This was put in 


2bXi - 


tZ 3 C i' V 
21 * iB . Z_ "^m-m 


mBJL 


by use of the relations 


*i " 


_ M 


2GAt 


B,- 


and 


Ci 


-2inC ± ' 


(97) 


From equations (20a), (22a), and (23a) and by use of the relations 
applicable to the circular tube of same chord and skin cross-section 
area. 
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From equations (22a) and (I 4 . 8 ) to (50), it may be shown that for 
the elliptic cross section 



(51) 

Similarly, from equations (20a), (3ii), (1+9 ) ^ and (50), the normal 
stress o^.' can be obtained as 


■ (t)(i)RH 


2 x: 


2 b 


sin 


2nns 


(52) 


Also, the shear stresses may be expressed as 



Airfoil-Shape Cross Section 

In order to determine the coefficients C n ' and C 0 ' of equa- 
tion (U), the contour of the airfoil box must first be expressed in the 
form r/l = f(s/l).. In this equation r is measured from the shear 
center of the box. The shear center of the box was determined by the 
assumption that the ct x stresses due to bending may be computed by the 
formula a x = Nyyi z where 

N bending moment 

y' distance from neutral axis to center of skin 

I z moment of inertia of skin about neutral axis (chord line for 

symmetrical wing) 
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Results for angles of twist 0/(M/lG) are plotted in figure 6(a) 
for rectangles with p = 0.05, 0.10, and 0.15. Values plotted were 
obtained using both values of obtained from equations (l8) and (19). 

For the rectangle p * 0.10, the results obtained by Kerman and Chien 
for the rectangle with p = 0.10 are plotted for comparison. The plots 
are of interest from several viewpoints. 

(1) The effect of end restraint on angle of twist disappears sub- 
stantially when x/2b is approximately 1.0, a distance of one chord 
length from the restrained end. (The chord here is considered the long 
side of the rectangle.) 

(2) Within small differences, the results obtained for the rectangle 
with p = 0.10 by Karman and Chien check with the results obtained 
using the value of X obtained from equation (18). 

(3) The series used for the computation of o/(M/l 0 G) of 0/(M/lG) 
converge with fair rapidity so that an excessive number of terms need 
not be taken. • 

(1+) Because the effect of the restraint is confined principally to 
the region near, x = 0, it is reasonable that the value of Xj_ satis- 
fying the basic equations in this region are the more logical ones to 
use. 

(5) The method of representing the contour of the box ty a Fourier 
series gives good results for the angles of twist of rectangular cross 
sections. 

(6) It should be noted here that the Karman-Chien values given are 
not necessarily exact, as these were also obtained by the use of a finite 
number of terms of an infinite series. 

In figure 6(b), results for angles of twist ©/(M/lG) are plotted 
for rectangles with P = 0.05, 0.10, 0.15, 0.20, 0.30, and O.Jj.0. Results 
for angles of twist 0y/(M/l o G) for rectangles with the same p values 
are shown in figure 7. The values plotted in both of these figures were 
based on X Q . 

• 

The calculations for a x ' gave values which are listed in table 8. 
These are plotted in figures 8 to 13. The series for a x ' given in 
equation (35) may be shown to be divergent at the corners at x = 0. 

This indicates an infinite stress (or the existence of a concentrated 
reaction) at the restrained end at the comers. Furthermore,' the series 
converges less rapidly at the comers of the box than elsewhere. This 
result does not agree with the results shown in figure 5 of reference 1, 
in which Karman and Chien show finite max im um stresses at the comers of 
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rectangular boxes. The series used by Karman and Chien may also be 
shown to be divergent at the comers at the restrained end. The result 
obtained by them was evidently obtained by use of only a finite number 
of terms in their series. In figures lit and 15> are shown comparisons 
of values of ct x ' as obtained in this report with values of ct x ' 
obtained by Karman and Chien and plotted in figure It of reference 1. 

The series for the shear stress (equation (3?)) is slowly conver- 
gent, particularly near the comers of the rectangles. However, ax. the 
restrained end the shear stresses are easily computed. At the restrained 
end x/2b = 0, 


T 0 1 - p(i)( m/i 0 g) 


( 61 ) 


For the rectangle, since r has only two values, this may be 
simplified to give the following: 

When r = b and x/2b = 0 


t _ 1 - ft 
T o " P 

When r = a and x/2b = 0 


(6la) 


At the unrestrained end 
stress is given by 


In table 9 are shown the results obtained from calculations using 
the first I4.O terms of the series. Because of the slow convergence the 
values are irregular. However, since the maximum shear values at the 
restrained end and the minimum values at the free end may be determined 
exactly, this difficulty in obtaining accurate values of shear stresses 
at intermediate stations was not felt to justify the extension of the 
computations to a larger number of terms. 

It appears from the study of rectangular sections that the existence 
of sharp comers is a major cause of high normal and shear stresses. 

This suggests that the possibility of rounding comers is worthy of 
investigation. 


t _ 


1.0 


(6lb) 


x/2b = 00 it may be shown that the shear 


T_ _ 1 

^O 


2ft 


(62) 



MCA TN 2600 


23 


Boxes of Elliptical Cross Sections 

Determination of coefficients .- In order to determine the basic 
series coefficients B n and C n , it is necessary to evaluate the 

integrals of equations (1;5) and (If?) . This was done by two methods. 

The first consisted of assuming values of 0 (equation (1*2)), computing 
values of (s /Z) and dr/ds from equations ( 38 a ), (lila), and (i+U3 and 
then evaluating the integrals of equations (1;5) and (I 4 . 7 ) by a numerical 
integration using a trapezoidal mile. The second method of evaluating 
these integrals consisted of plotting the functions (dr/ds) sin (2nits/i ) 

and (s/Z - 0/2n ) sin (2nns/l) against s/Z and integrating graphi- 
cally by the use of a planimeter. In table 10 are listed values obtained 
for the coefficients by these two methods for an ellipse with an a/b 
ratio of O.O 87 I. 

Computation for angle.s of twist and stresses.- In table 11 are 
given value s of 9 / (M/l 0 G) for this ellipse. Columns (l) and (2) were 
obtained by using coefficients C n and B n and obtained by calcu- 

lations. Values given in column (l) were obtained by using the first 
10 coefficients. Values given in column (2) were obtained by using the 
first 20 coefficients-. Values in column (3) were obtained by using the 
first 10 coefficients obtained by the graphical determination. 

At the unrestrained end, the simple torsion theory can be applied^ 
This leads to a value of 



For the ellipse with a/b = O.O 87 I, this equation gives 

= 55-00. A comparison of this value with the values in table 11 

shows that a 20-term series, whose coefficients were obtained by calcu- 
lation, gives a very good approximation at infinity. 

For the ellipse with a/b = 0.1738, calculations were based on a 
20-term series using coefficients obtained by the approximate computa- 
tional method of determining the integrals in equations (lj.5) and (I 4 . 7 ) . 

For all other ellipses calculations were based upon a 10-term series. 

In tables 12 to 1 5 are given coefficients C n , Bn, and (2\ n b) 0 

that were used in the computation of 0/(M/l o G), a x ' , and t/t 0 for 
the ellipses investigated. 
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In table 16 values are given for angles of twist for the various 
ellipses. These are plotted in figures 16 and 17. Both 9 /(m/i o G) and 

0/(H/lG) are plotted. 

Values of a x ' were computed using the coefficients given in 
columns (2), (3), and (!;) of table 10 and those given in tables 12 to 
13. These c^' values are given in table 17 and are plotted in fig- 
ures 18 to 22. 

Values of t/t q are given in table 18. 


Box of Airfoil Cross Section 


The airfoil was drawn to a large scale. The coordinate system 
shown in figure 3 was adopted. The shear center was determined to be 
36.8 units from the leading edge, the chord being taken as 63.0. The 
length of the periphery was obtained by summing up small distances As 


along the contour. This length l = ) As = llj.2.1. The enclosed area 

r 

of the cross section was found to be 663.0. -The value of i/t = '€) ds 


p v 

was obtained by a summation of (r^As) over the periphery. This gave 
I/t = 17,270, In table 19 are given values of r/l against s/l 
obtained from measurements on the airfoil. 


Coefficients C 0 ' and C n ' were obtained from equations (3U) and 
(33). The integrals were evaluated by replacing the integrals by sum- 
mations. In table 20 are given the coefficients obtained. Twenty coef- 
ficients were used in the calculations. Values of r/l were then com- 
puted at various values of s/l. These are listed in table 19 for 
comparison with the measured values. These two sets of r/l against 
s/l values from table 19 are plotted in figure 23. The accuracy of the 
series is fairly good except near, the points of discontinuity of the 
profile . 

The coefficients B n were determined by equation (36), the inte- 
grals being replaced by summations. Values of these coefficients are 
listed in table 20. 

Coefficients 2X n b were computed from equation (37). These are 
also listed in table 20. 

In table 21 and figure 2k are given the results of computation for 
angles of twist. These were computed from equation (38). Angles of 
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twist rapidly approach a constant value as the parameter x/2b is 
increased. The effect of the restraint at the inner end of the box on 
angles of twist disappears almost entirely in one chord length from the 
restrained end. 

The results of the calculations for the normal stresses <r x ' are 
given in table 22 and figure 25. The infinite stress at the comers of 
the cross section at the restrained end is again shown. In a distance 
of one chord length from the restrained end, however, the effect of the 
restraint has been reduced so that the normal stresses are negligibly 
small. The series used in computing ct x ! , converges with fair rapidity 
so that use of 20 coefficients gives good accuracy, except right at the 
restrained end. Here the series converges less rapidly and is divergent 
at the sharp comers. 

Table 23 gives the computed values of t/t 0 for the airfoil sec- 
tion. The series giving these values does not converge rapidly near the 
restrained end. The values given in this* table are based on 20 coef- 
ficients in the series. This is insufficient for good accuracy. 


CONCLUSIONS 


An approximate method for solving the basic equations of Karman and 
Chien was determined and applied to the calculation of angles of twist 
and stresses in torsion boxes of rectangular, elliptical, and airfoil 
cross section. From the results of the application of this method, the 
following conclusions may be drawn: 

1. The methods given here for the approximate solution of the 
Karman-Chien equations are applicable to torsion boxes having one degree 
of symmetry of the cross section. The series used to obtain Angles of 
twist converge fairly rapidly. The results check well for the rectan- 

. gular cross sections with those given by Karman and Chien. The series 
for ct x stresses converges less rapidly, particularly at the restrained 

end. At the comers of the restrained end, the series diverges, indica- 
cating the necessity for concentrated loads at these points. 

2. The theory gives results which indicate tendencies that might 
well be considered in design. The infinite stresses at comers at the 
restrained end suggest the possibility of high stress concentrations at 
these points. It would be expected that, even with zero circumferential 
strains, the computed values for angles of twist and stresses would be 
modified appreciably by a local yielding in the neighborhood of the 
restrained ends, particularly at sharp comers. 
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3. The studies also indicate that a wing section composed of 
approximately square boxes might be of some advantage in reducing angles 
of twist and decreasing normal stresses due to torsion. 


Oregon State College 

Corvallis, Ore., March 28, 19f>0 
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T ABLE X.- SERIES COEFFICIENTS FOR RECTANGLE WITH 0 = 0.05 


n 

C n 

Bn 

( 2 V>). 

( 2 x n b )o 

2 

-0.2813 

-0.07507 

5.71 

2.88 

1* 

-.5559 

-.03708 

11.1*3 

5.75 

6 

-.8169 

-.021*21 

17.17 

8.62 

8 

-1.0575 

-.01763 

22.96 

11.50 

10 

-1.2723 

-.01358 

28.78 

H*. 37 

12 

-1.1*559 

-.01078 

31*. 65 

17.25 

ll* 

-1.6028 

-.00873 

1*0.58 

20.12 

16 

-1.7112 

-.00713 

1*6.55 

23.00 

18 

-1.7770 

-.00586 

52.57 

25.87 

20 

-1.7998 

-.001*80 

58.62 

28.71* 

22 

-1.7767 

-.00392 

61*. 68 

31.62 

2 b ' 

-1.7122 

-.00317 

70.81 

31*. 1*9 

26 

-1.6031* 

-.00253 

76.93 

37.37 

28 

-1.1*560 

-.00198 

83.01* 

1*0.21* 

30 

-1.3263 

-.00151 

89.11 

1*3.12 

32 

-1.0570 

-.00111 

9^.23 

1*5.99 

31* 

-.8178 

-.00075 

101.29 

1*8.87 

36 

-.5562 

-.0001*6 

107.33 

51.71* 

38 

-.2816 

-.00021 

113.35 

51*. 61 

bO 

0 

0 

119.32 

57-1*9 

b 2 

.2822 

.00017 

. 125.29 

60.36 

bb 

.3514; 

.00031 

131.20 

63.21* 

b 6 

.8175 

.0001*1 

137.12 

66.11 

1*8 

1.0581 

.0001*9 

11*3.03 

68.98 

50 

1.2723 

.00051* 

11*8.91* 

71.86 

52 

1.1*565 

.00057 

151*. 85 

71*. 71* 

51* 

1.601*0 

.00059 

160.78 

77.61 

56 

1.7127 

.00058 

166.71 

80.1*8 

58 

1.7775 

.00056 

172.67 

83.36 

60 

1.8011 

.00053 

178.6I* 

86.23 

62 

1.7755 

.0001*9 

181*. 61 

89.ll 

6 b 

1.7122 

.0001*5 

190.63 

91.98 

66 

1.5653 

.00039 

196.62 

9l*. 86 

68 

1.1*31*8 

.00031* 

202.61* 

97.73 

70 

1.2701* 

.00028 

208.66 

100.60 

72 

1.0581 

.00022 

211*. 67 

103.1*8 

lb 

.8179 

.00016 

220.68 

106.35 

76' 

.7875 

.00010 

226.61* 

109.23 

78 

.2792 

.00005 

232.67 

112.10 

80 

0 

0 

233.61* 

111*. 98 
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TABLE 2.-- SERIES COEFFICIENTS FOR RECTANGLE WITH p = 0.10 


n 

C n 

■ ^ 

(X b L 

(V). 

2 

-0.1+91+1 

-0.06957 

5.25 

3.1+7 

1+ 

-.9397 

-.03309 

10.57 

6.91+ 

6 

-1.2936 

-.02021+ 

16.05 

10.1+1 

8 

-1.5208 

-.01338 

21.67 

13.88 

10 

-1.5989 

-.00906 

27.1+3 

17.35 

12 

-1.5208 

-.00591+ 

33.29 

20.81 

H+ 

-1.291+2 

-.00371 

39.18 

21+.28 

16 

-.9395 

-.00207 

1+5.01+ 

27.76 

18 

-.1+91+0 

-.00086 

50.85 

31.21+ 

20 

0 

0 

56.50 

31+.68 

22 

.1+91+6 

.00057 

62.17 

38.18 

21+ 

.91+05 

.00092 

67.68 

1+1.62 

26 

1.2932 

.00108 

73.29 

1+5.11 

28 

1.5210 

.00109 

78.89 

1+8.58 

30 

1-5976 

.00100 

81+.53 

52.05 

32 

1.5213 

.00081+ 

90.21 

55.52 

3 k 

1.2939 

.00063 

95.81+ 

58.92 

36 

• 91+05 

.0001+1 

101.72 

62.1+8 

38 

.1+91+0 

.00019 

107.28 

65.86 

1+0 

0 

0 

113.00 

69.36 

1+2 

-.1+932 

-.00016 

118.78 

72.92 

kk 

-.5995 

-.00028 

121+.28 

76.32 

1+6 

-1.2913 

-.00035 

129.88 

79.79 

1+8 

-1.5193 

-.00037 

135. 61+ 

83.35 

50 

-1.5991 

-.00036 

H+1.13 

86.73 

52 

-1.5218 

-.00031 

H+6.85 

90.23 

51+ 

-1.2920 

-.00025 

152.59 

93.73 

56 

-.9390 

-.00018 

158.25 

97.16 

58 

-.1+951+ 

-.00008 

I6I+.00 

100.66 

60 

0 

0 

169.70 

101+.16 

62 

.1+91+5 

.00007 

175.20 

107.51+ 

61+ 

.91+05 

.00013 

180.95 

111. 10 

66 

1.2932 

.00017 

186.1+1 

1H+.1+7 

68 

1.1+979 

.00018 

- 191.88 

U7.85 

70 

1.6001 

.00018 

197.68 

121.1+1 

72 

1.5238 

.00016 

203.1+0 

121+.91 

71+ 

1.2919 

.00013 

208.92 

128.28 

76 

.91+02 

.00009 

211+.96 

131.97 

78 

.1+91+7 

.00005 

221.00 

135.65 

80 

0 

0 

226.00 

138.72 
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TABLE 3.- SERIES COEFFICIENTS FOR RECTANGLE WITH p « 0.15 


n 

c n 

/ B n 

(2V)„ 

( 2 V>)o 

2 

-0.63551* 

-0.06320 

U.89 

3.85 

1* 

-1.13291 

-.02815 

9.97 

7.71 

6 

-1.3821*8 

-.01527 

15.32 

11.56 

8 

-1.33086 

-.00827 

20.88 

15.1*1 

10 

-.98973 

-.001*02 

26.1*8 

19.27 

12 

-.U3257 

-.00119 

31.96 

23.10 

u* 

.21892 

.0001*1* 

37.39 

27.00 

16 

.82293 

.00128 

1*2.61 

30.83 

18 

1.13266 

.00139 

U7.85 

31*. 65 

20 

1.39918 

.00139 

53.19 

38.55 

22 

1.21*758 

.00102 

58.59 

1*2. 1*0 

21* 

.8211*2 

.00057 

63.91* 

1*6.20 

26 

.21880 

.00013 

69.1*0 • 

50.10 

28 

— . Jj-3257 

-.00022 

71*. 71* 

53.96 

30 

-.98910 

-.0001*1* 

79.93 

57.75 

32 

-1.33236 

-.00053 

85.28 

61.65 

3 h 

-1.38361 

-.0001*7 

90.63 

65.51 

36 

-1.13266 

-.00035 

95.92 

69.30 

38 

— .63U78 

-.00018 

101.1*0 

73.22 

Uo 

0 

0 

106.90 

77.17 

k 2 

.63566 

.00015 

112.08 

80.92 

1*1* 

1.13015 

.00023 

117.31* 

81*. 75 

b 6 

1.38085 

.00026 

122.71 

88.65 

1*8 

1.32935 

.00023 

128.12 

92.55 

50 

.98960 

.00016 

133.26 

96.23 

52 

.1*3106 

.00006 

138.80 

100.20 

51* 

-.2201*3 

-.00003 

11*1*. 00 

103.95 

56 

-.82293 

-.00011 

11*9.1*7 

107.92 

58 

-1.21*570 

„ -.00015 

151*. 81* 

111.82 

60 

-1.39793 

-.00015 

160.02 

115.58 

62 

-1.13301* 

-.00011 

165.1*5 

119.1*7 

61* 

-.81992 

-.00008 

170.97 

123.1*1* 

66 

-.21967 

-.00002 

176.10' 

127.13 

68 

.1*21*97 

.00001* 

181.50 

131.02 

70 

.98910 

.00008 

186.86 

131*. 92 

72 

1.33387 

.oooio 

191.91* 

138.60 

7l* 

1.1*0810 

.00010 

197.31* 

11*2.50 

76 

1.13H5 

.00008 

202.96 

11*6.55 

78 

.63679 

.00001* 

208.00 

150.16 

80 ' 

0 

0 

213.50 

151*. 13 
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TABLE 4.- SERIES COEFFICIENTS FOR RECTANGLE WITH p - 0.20 


n 

C n 

% 

( 2 M 0 

2 

-0.70362 

-O.0339I 

4.04 

4 

-1.14143 

-.02261 

8.07 

6 

-1.H033 

-.OlOOli 

12.11 

8 

-.70387 

-.00349 

. 16.16 

10 

0 

0 

20.19 

12 

.70612 

.00133 

24.23 

Hi 

.1.14120 

.00184 ' 

28.27 

16 

1.14043 

.00142 

32.30 

18 

.70337 

.00069 

36.33 

20 

0 

0 

40.38 

22 

-.70600 

-.00046 

44.42 

2Li 

-l. 14093 

-.00063 

48.46 

26 

-1.11033 

-.00033 

32.30 

28 

-.70160 

-.00029 

36.31 

30 

0 

0 

60.36 

32 

.70337 

.00022 

64.61 

34 

1.14020 

.00031 

68.63 

36 

1.UO70 

.00028 

72.69 

■a .• 

.70637 

' .00013 

76.73 


0 

0 

80.77 


-.70688 

-.00013 

84.80 


-1.HO20 

-.00019 

88.84 

46 

-1. 14108 

-.00017 

92.88 

48 

-.70337 

-.00010 

96.92 

3o 

0 

0 

100.96 

32 

.70337 

.00008 

104.99 

34 

1.14283 

.00012 

109.04 

36 

1.13944 

.00011 

113.07 

38 

.70298 

.00006 

117.11 

60 

0 

0 

121.13 

62 

- . 7014.7U 

-.00006 

123.18 

6U 

-l.llilli3 

-.00009 

129.23 

66 

-1.14396 

-.00003 

133.26 

68 

-.69U23 

-.00003 

137.31 

70 

0 

0 

141.34 

72 

.70337 

.00004 

143.38 

74 

1.13836 

.00007 

149.42 

76 

1. Ili070 

.00007 

133.43 

78 . 

.70337 

.00004 

137.30 

80 

0 

0 

161.33 
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TABLE SERIES COEFFICIENTS FOR RECTANGLE WITH p = 0.30 


n 

C n 

B n 

0V>)o 

2 

-0.61i722 

-O.O3907 

l *. oi * 

k 

-.76088 

-.0111*9 

8.07 

6 

■ -.21*718 

-.00166 

12.10 

8 

.1*7075 

.00178 

16 . ll * 

10 

.80007 

.00193 

20.17 

12 

■ .1*7025 

.00079 

21*. 21 

Hi 

-.21*706 

. -.00031 

28.23 

16 

-. 7616 k 

-.00072 

32.27 

18 

-.61772 

-.0001*9 

36.30 

20 

0 

0 

1*0.31* 

22 

.61*659 

.00032 

1*1*. 38 

2 ii 

-.76261* 

.00032 

1 * 8 . 1*0 

26 

-.21*655 

.00009 

52 . 1 * 1 * 

28 

, -.1*7125 

-.00015 

56.1*7 

30 

-.79882 

-.00021 

60.51 

32 

-.1*7025 

-.00011 

61*. 51* 

3 k 

.21*768 

.00005 

68.57 

36 

.76189 

. 00011 * 

72.61 

38 

.61*910 

.00011 

76 . 61 * 

liO 

0 

0 

80.68 

1*2 

-.61*885 

-.00009 

81*. 71 

iili 

-. 76261 * 

-.00010 

88.71* 

li 6 

-.21*81*1* 

-.00003 

92.78 

1*8 

.1*7025 

.00005 

96.81 

30 

.80111 

.00008 

100.81* 

32 

.1*7025 

. 00001 * 

10 l *.88 

3 k 

-.21*756 

-.00002 

108.91 

36 

-75963 

-.00006 

112.95 

38 

-.61*835 

-.00005 

116.98 

6 o 

0 

0 

121.01 

62 

.65023 

. 00001 * 

125.05 

6 k 

.75963 

.00005 

129.08 

66 

.21*869 

.00001 

133.12 

68 

-. 1 * 6281 * 

-.00002 

137. ll* 

70 

-.80007 

-. 00001 * 

Ha. 18 

72 

-.1*7025 

-.00002 ' 

11*5.22 

7 U 

. 21*630 

.00001 

H*9.25 

76 

.76365 

.00003 

153.29 

78 

.61*659 

.00002 

157.31 

80 

0 

0 

161.35 
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TABLE 6 „- SERIES COEFFICIENTS FOR RECTANGLE WITH 0 = 0 . 4 ° 


n 

c n 

B n 

( 2 V>)o 

2 

- 0.38044 

- 0.02007 


4 

-.23412 

-.00310 


6 

.23512 

.OOI38 


8 

.38032 

.00126 

14.88 

10 

0 

0 

18.61 

12 

-.38057 

-.00056 

22.33 

li+ 

-.23563 

-.00025 

26.15 

16 

.23512 

.00019 

29.11 

18 

.38094 

.00025 

33.4 9 

20 

0 

0 

37.21 

22 

-.37994 

-.00016 

40.94 

24 

-.23512 

-.00009 

44 066 

26 

.23512 

.00007 

48.37 

28 

.38157 

.00010 

52.10 

30 

0 

0 

55.82 

32 

-.■37981 

-.00008 

59.54 

3 k 

-.23487 

-.00004 

63.27 

36 

.23512 

.00004 

66.98 

38 

.38182 

.00006 

70.70 

ko 

0 

0 

74.43 

h 2 

-.38245 

-.00005 

78.15 

h k 

-.23487 

-.00002 

81.87 

h 6 

.23399 

.00002 

85.59 

48 

137981 

.00004 

89.31 

50 

0 

0 

93.03 

52 

-.37881 

-.00003 

96.76 

54 

-.23399 

-.00002 

ioo „48 

56 

.23563 

j .00002 

104.19 

58 

.38245 

.00003 

107.92 

60 

0 

0 

111.64 

62 

-.38157 

-.00002 

115.36 

6 h 

-.23311 

-.00001 

119.08 

66 

' .23211 

.00001 

122.80 

68 

.37448 

.00002 

126.52 

70 

0 

0 

130.24 

72 

-.37981 

-.00002 

133.97 

74 

-.23701 

-.00001 

137.69 

76 

.23387 

.00001 

141.40 

78 

.38208 

.00001 

145.13 

80 

0 

0 

148.85 
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TABLE 7.- ANGLES OF TWIST 0/(M/l o G) AND 0/(M/lG) 

FOR RECTANGLES 

(a) For rectangles with p = 0.05, 0.10, and 0.15* 


X 

2 b 

Based on X^ 

Based on X Q 

9 

9 

0 

0 

I^G 

m/ig 

M/I 0 G 

ll/l'G 

p = 0.05 

0 

19.00 

1.00 

19.00 

1.00 

1/8 

81.08 

u .26 

68.11 

3.58 

1/h 

89.50 

U .71 

81.15 

u .27 

1/2 

91;. 6 0 

U .98 

8 9-59 

U .72 

3/k 

95.03 

5.oi 

92.1*8 

1; . 86 

1.0 

95 .1*2 

5.02 

93.97 

U .9 U 

2.0 

95.70 

5.01; 

95 A 2 

5.02 

p = 0.10 

0 

9.00 

1.00 

9.00 

1.00 

1/8 

19.25 

2.16 

17.07 

1.90 

1A 

21.91; 

2.1*3 

20.37 

2.26 

1/2 

23.61; 

2.63 

22.73 

2.53 

3A 

21;. 17 

■2.69 

23.62 

2.62 

1.0 

21;. 37 

2.71 

21;. 03 

2.67 

2.0 



21;. 53 

2.73 

2 U . U 8 

2.72 

p = 0.15 

0 

5.67 

1.00 

5-67 

1.00 

1/8 

8.79 

1.55 

8.23 

1 A 5 

iA 

9.68 

1.71 

9.30 

1.61; 

1/2 

10.1;8 

1.85 

10. 2 lj . 

1.81 

3A 

10.76 

1.90 

10.62 

1.87 

1„0 

10.89 

1.92 

10.80 

1.91 

2.0 

10.98 

1.9 U 

10.98 

1.9 U 
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TABLE 7.- ANGLES OF TILE ST 0/(m/I o g) AND 0/(M/lG) 
FOR RECTANGLES' - Concluded 
(b) For rectangles with p = 0.20, 0.30, and 0.1*0. 


X 

2b 

Based on \ 0 

0 

9 

M/I 0 G 

m/ig 

p = 0.20 

0 

1*.00 

1.00 

1/8 

U.98 - 

1.23 

lA 

5.1*2 

1.35 

1/2 

5.85 

1.1*6 

3/1* 

6.03 

1.51 

1.0 

6.12 

1.53 

2.0 

6.21 

1.55 


p =0.30 

. 

0 

2.33 

1.00 

1/8 

2.50 

1.07 

1 /k 

• 2.59 

1.11 

1/2 

2.67 

1.15 

3/1* , 

2.72 

1.17 

1.0 

2.71* 

1.18 

2.0 

2.77 

1.19 

p = 0.1*0 
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TABLE 8.- NORMAL STRESSES ct x ' FOR RECTANGLES 
BASED ON (2X n b) 0 


(a) p = 0.05. 
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TABLE 8.- NORMAL STRESSES a x ' FOR RECTANGLES 
BASED ON (2\ n b) 0 - Concluded 


(d) 0 = 0.20. 


X 

2b 

CM 

O 

0 

O 

II 

CQ|r~» 

7 = OoOli 

Is 

f ■ °-° 5 

n 

7 = 0.10 

if 

- = 0.20 

If 

0 

-o.Bk 

-1.18 



-1.28 

-0.58 

-0.12 

1/8 

-.38 

-.71; 

-0.81; 

-.81; 

-.53 

-.13 

lA 

-.23 

-.la 

— -1|7 

-.50 

-,1|2 

-.12 

1/2 

-.09 

-.16 

-.19 

-.21 

-.23 

-.09 

3A 

-.oil 

-.07 

— 

-.11 

-.12 

-.06 

1.0 

-.02 

-oOjU 

— 

-.OB 


-.01; 

2.0 

0 

0 

— 

0 


0 




(e) 0 = 

0.30. 



X 

2b 

7 = O.O3 

is 


f = 0.075 

m 


f - °- 20 

0 

-0.18 

-0.1|9 

. 

-0.1|0 

-0.18 


1/8 

-.16 

-.3k 

-O.liO 

-.33 

-.17 

■Kfl 

lA 

-.12 

-.22 

-.23 

-.2k 

-.15 

-.07 

1/2 

-.06 

-.10 

-.11 

-.12 

-.10 

-. OB 

3A 

-.03 

-.05 

-.06 

-.06 

-.06 

-.03 

1.0 

-.01 

-.03 

— 

-.03 

-.03 

-.02 

2.0 

0 

0 

— 

0 

0 

0 

(f) p = 0.1|0. 

X 

mm 

3 = 0.06 


I 


£ = 0.20 

2b 

BBlfil 

l 

* fgf 

mmmm 

wmmm 

1 

0 

-0.05 

1 


-0.16 

-0.09 


1/8 

-.oh 


-0.16 

-.A 

-.09 

-.01; 

lA 

-.Oli 

-.07 

-.10 

-.10 

-007 

-.03 

1/2 

-.02 

-.0l| 

-.05 

-.OB 

-.05 

-.02 

3A 

-.01 

-.02 

-.03 

-.03 

-.03 

-.02 

1.0 

-.01 

-.01 


-.02 

-.02 

-.01 

2.0 

0 

0 

— 

0 

■ 

0 

0 
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TABLE 9o- SHEAR STRESSES t/t q FOR RECTANGLES 
BASED ON (2A n b) 0 


(a) p = 0.0 5. 


1=0 

l 

f = 0.003 

L 

<r-J |C/) 

II 

o 

o 

- = 0.02 
1/ 

|co 

II 

0 

• 

H 

O 

|co 

II 

0 

0 

rv> 

0 

f ■ °- 25 

19.00 

21.11 

17.21; 

15.68 

15.27 

15.23 

1U.02 

19.00 

1.03 

-1.61 

-2.85 

-3.22 

-3.25 

-3.56 

19.00 
19.18 
19.92 
19. hO 
19.15 
19.15 
18.87 

1.00 

29.32 

26.31; 

20.77 

18.19 

16.89 

15.30 

1.00 

il.79 

6.72 

8.70 

9.U5 

9.68 

9.61 

1.00 

1;.08 

5.28 

6.71 

7.65 

8.31 

9.l£ 

1.00 

U.05 

5.26 

6.59 

7.51 

8.19 

9.1;2 


(b) p = 0.10. 


X 

2b 

s = 0 

l 

f = 0.01 

If 

- = 0.02 

f = °-° 5 

L 

- = 0.10 
l 

- = 0.20 
l 

0 

9.00 

9. 00 ■ 

9.00 

1.00 

1.00 

1.00 

1/8 

8.95 

7.58 

3.23 

2.87 

2.35 

2.23 

1/h 

7.18 

6.32 

3.29 

5.75 

3.6 9 

2.77 

1/2 

6.19 

5.52 

3.08 

5.66 

1;.60 

3.63 

3/h 

5.86 

5.26 

2.92 

5-25 

4.87 

1;. Ql; 

1.0 

5.75 

5.15 

2.81; 

U.97 

h-93 

h.3$ 

2.0 

5.62 

5.03 

2.7U 

l;.6l 

4.89 

1*.77 






(c) 

p = 0.15. 




JL 

2b 

- = 0 
l 

- = 0.01 
l 

^ = 0.02 
l 

- = 0.05 

l 

- = 0.10 
l 

- = 0.20 
l 


0 

5.67 

5.67 

5.67 

1.00 

1.00 

1.00 

1.00 

1/8 

5.39 

5.07 

5A9 

3.50 

1.98 

1.67 

1.65 

lA 

u.ia 

1; . 28 

U.9i 

3.79 

2.60 

2.07 

2.01; 

1/2 

3.87 

3.67 

UA7 

3.52 

3.12 

' 2.58 

2.5U 

3A 

3.65 

3A6 

h.3l 

3.27 

3.25 

2.89 

2.86 

1.0 

3.55 

3.37 

h.23 

3.12 

3.26 

3.07 

3.05 

2.0 

3A7 

3.29 

U.15 

2.95 

3.60 

3.29 

3.20 
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TABLE 9.- SHEAR STRESSES x / t q FOR RECTANGLES 
BASED ON (2)^ ) Q - Continued 


(d) p = 0.20. 



^NACA^ 
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TABLE 9.- SHEAR STRESSES tt/t 0 FOR RECTANGLES 
BASED ON (^^b) 0 - Concluded 


(f) p = 0.U0. 


X 

2b 

<x |oa 
■II. 
O 

£ = 0.02 
z 

f = 0.0625 

V 

- = 0.10 
z 

^21 

- = 0.20 
z 

f=0. 

0 

1.50 

1.50 

1.50 

1.50 

1.00 

1.00 

1.00 

1.00 

1/8 

1.24.7 

1.24-6 

1.U3 

1.51 

1.00 

1.05 

l.ob 

1.0b 

lA 

1.1*2 

1.1a 

1.37 

1.52 

1.00 

1.10 

1.09 

1.09 

1/2 

1.36 

1.30 

1.33 

1.52 

1.00 

1.17 

1.15 

1.15 

3A 

1.31 

1.30 

. 1.29 

1.51 

.99 

1.20 

1.19 

1.18 

1.0 

1.29 

1.29 

1.28 

1.52 

1.00 

1.23 

1.22 

1.22 

2.0 

1.26 

1.25 

1.26 

1.51 

.99 

1:25 

1.2b 

1 1,25 
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TABLE 10.- SERIES COEFFICIENTS FOR ELLIPSE 
WITH a/b = 0.0871 


(1) 

(2) 

(3) 

(U) 

(5) 

(6) 

(7) 


Coefficients by 

Coefficients by 

n 

calculation 


graph integration 


C n 

V 

( 2 vOo 

C n 

Bn 

( 2X n b )o 

2 

0.2501 

0. 01*677 

i*.9l* 

0.251 

0.01*71 

5.12 

U 

-.3625 

-.01639 

10. ll* 

-.381* 

-.0168 

10.05 

6 

.1*270 

.0081*7 

11*. 71 

.1*1*1* 

.0086 

15.05 

8 

— .I4.876 

-.00550 

19.51* ■ 

-.512 

-.0063 

20.67 

10 

.5338 

.00391* 

21*. 52 ' 

.582 

.0036 

21*. 23 

12 

-.5725 

-.00290 

29.1*1 

-.595 

-.0031* 

31.31 

H* 

.6032 

.00226 

31*. 50 

.678 

.0020 

33.20 

16 

-.6350 

-.00181 

39.18 

-..611* 

-.0011 

35.80 

18 

.5868 

.0011*7 

1*1*. 10 

.672 

.0010 

1*2.85 

20 

-.6821 

-.00116 

1*8.55 

-.71*9 

-.0012 

1*8.91* 

22 

.6992 

.00105 

53.80 

.768 

.0015 

— 

2l* 

-.1192 

-.00092 

58.98 

-.791* 

-.0010 

— 

26 

.7359 

.00081 

61*. 07 

.791* 

.0009 


28 

-.71+77 

-.00069 

68.31 

-.838 

-.0008 

— 

30 

.7601 

.00058- 

71*. 19 

.813 

.0007 


32 

-.7669 

-.0001*8 

78.91* 

-.851 

-.0006 

— 

3 k 

-771+5 

.00052 

83.87 

.838 

.0005 

— 

36 

-.7831 

-.0001*7 

87.96 

-.787 

-.0001* 

— 

38 

.7882 

.00016 

93.99 

.870 

.0001* 

— 

1+0 

-.7957 

-.00105 

97.78 

-.857 

-.0005 

— 







TABLE 11.- ANGLES OF TWIST e/(M/l 0 G) AND e/(M/lG) FOR ELLIPSE 

WITH a/b = 0.0871 



( 1 ) 

( 2 ) 

(3) 

z 

2b 

Obtained, by use of 
first 10 coefficients 
of table 10, columns 
(2) , (3), and A) 

Obtained by use of 
first 20 coefficients 
of table 10, columns 
(2), (3), and (1*) 

Obtained by use of 
first 10 coefficients 
of table 10, columns 
( 5 ), ' ( 6 ), and ( 7 ) 


8 


0 

0 

0 

§ 


K/IqG 

m/ig 

M/I 0 G 

m/ig 

M/I 0 G 

M/IG 

0 

3^-90 

1.00 

31*. 90 

1.00 

31 *. 90 

1.00 

1/8 

1*5.02 

1.29 

1*7-29 

1.36 

1*1 *. 67 

1.28 

1/li 

48.16 

1.38 

50.51* 

■1.1*5 

1*8.16 

1.38 

1/2 

30.50 

1.1*5 

52.98 

1.52 

50.61 

1.1*5 

3A 

51.23 

1-47 

53.89 

i.51* 

51:1*8 

1.1*7 ■ 

1.0 

51.58 

1.1*8 

54.27 

1.56 

51.86 

1.1*9 

2.0 

5?. 11 

1.1*9 

54.58 

1.56 

52.11* 

1.1*9 

00 

52.15 

1.1*9 

51 *. 62 

1.56 

52.18 

1.1*9 


asssaga 
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TABLE 12.- SERIES COEFFICIENTS FOR 
ELLIPSE WITH a/b = 0.1738 


NACA TN 2600 


n 


Bn 

( 2 vOo 

2 

0.1*378 

o.ol*l*o5 

3.22 

k 

-.3791 

-.011*37 

10.39 

6 

.631*8 

.00733 

13-36 

8 

-.6892 

-.ool*33 

20.76 

10 

.7223 

.00291 

23.86 

12 

-.7227 

-.00201 

31.38 

Hi 

.7122 

.0011*3 

36.1*0 

16 

-.6982 

-.00109 

39.33 

18 

.681*8 

.00083 

1*3-38 

20 

-.661*3 

-.00068 

1*9.38 

22 

.614*3 

,00032 

37.1*1* 

2k 

-.6198 

-.oool*l* 

63-09 

26 

.39U; 

.00036 

68.1*6 

28 

-.3679 

-.00028 

72.89 

30 

.31*20 

.00021; 

78.31* 

32 

-.331*9 

-.00020 

82.89 

3h 

.1*838 

.00016 

87.93 

36 

-.1*382 

-.00012 

91.31* 

38 

.1*31*6 

.00012 

99.71 

1*0 

-.1*091 

-.00008 

100.13 
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TABLE 13.- SERIES COEFFICIENTS FOR 
ELLIPSE WITH a/b = 0.2588 ' 


n 

C n 

Bn 

( 2x nt) 0 

2 

o.5itOit 

0.0388 

5.31 

it 

-.6i|39 

-.Ollit 

10.60 

6 

.6560 

.0056 

16.02 

8 

-.6317 

-.0028 

21.22 

10 

.5873 

.0017 

26.58 

12 

-.5Uoi 

-.0011 

31.90 

lit 

.it868 

.0008 

37.52 

16 

-.it333 

-.OOQit 

itl.60 

18 

.3663 

.0003 

it7.it8 

20 

-.36it9 

-.0002 

51.72 











NACA TN 2600 


k6 


TABLE lit.- SERIES COEFFICIENTS FOR 
ELLI PSE WITH a/b = 0.5000 


n 

c n 

Bn 

K b ) 0 

2 

0.5287 

0.0261t0 

5.0it 

It 

-.3989 

-.OOlt98 

10.15 

6 

.2753 

.0015U 

15.10 

8 

-.17U8 

r.ooo^lt 

20.09 

10 

.0982 

.00018 

25.22 

12 

-.061t3 

-.00009 

30.33 

lit 

.0386 

.OOOOlt 

35.1a 

16 

-.0231 

-.00002 

ltO.60 

18 

.0136 

.00001 

U6.60 

20 

-.0070 

0 

51.85 


TABLE 15.- SERIES COEFFICIENTS FOR 
ELLIPSE WITH a/b - 0.8660 
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TABLE 16.- ANGLES OF TWIST b/(K/I q g) AND ©/(M/lG) FOR 
ALL ELLIPSES FROM (21 n b) Q 



0 3k.90 

1/8 1+7-2 9 
i/h 50.5U 
1/2 52.98 
3 A 53.89 

1.0 51*. 27 

2.0 51**58 

00 51**62 


1.00 

1.36 

1.1*5 

1.52 


,1738 I 


0 


M/IG 


10.97 1.00 

12.88 1.17 


13.59 
11*. 16 


1.51* 11*. 38 

1.56 ll*.l*8 
1.56 11*. 55 
1.56 H*. 55 


1.21* 

1.29 

1.31 

1.32 

1.33 
1.33 




.2588 £ = 0.5000 £ 
D D 


0 9 0 0 9 


M/lG M/I 0 G M/IG . M/I 0 G M/IG 


1.00 

1.10 

l.ll* 

1.18 

1.20 

1.20 

1.21 


1.00 

l.Ol* 

1.07 

1.09 

1.11 

1.11 

1.11 

1.11 








































NACA TN 2600 


1*8 


TABLE 17.- NORMAL STRESSES ct x ' FOR ALL 
ELLIPSES FROM (2A n b) 0 


(a) a/b = 0.0871. 


X 

2b 

s 1 

T ~ 16 

s 1 

T = S 

s _ 3 
l 16 

s _ 7 
I 32 

0 

0.1*806 

1.0609 

2.33914. 

3.3680 

1/8 

.1*1*1*!* 

.9727 

1.3020 

1.3068 

lA 

.3778 

.71M 

.9179 

.6832 

1/2 

.2263 

.3872 

.31*37 

.208k 

3A 

.1229 

.1913 

.1303 

. O8I4.2 

1.0 

.062I* 

.0936 

.0699 

.0383 

2.0 

.0037 

.0032 

.0037- 

.0020 


(b) a/b = 0.1738. 


X 

2b 

s _ 1 

, T ~ 16 . 

. s _ i 
l 8 

^ | cn 
II 

s = J_ 
l 32 

0 


0.6261 

1.1263 

1.3797 

1/8 


.1*833 

.71*68 

.681*1* 

lA 

mmmEna' ' 

.3728 

.1*363 

.3069 

1/2 

.1108 

.1830 

.1632 

.0933 

3A 

.036U 

.0863 

.0673 

.0377 

1.0 

.0278 

.01*09 

.0302 

.0163 

2.0 

.0012 

.0016 

.0012 

.0008 


(c) a/b = 0.2388. 


X 

2b 

s = 1 

I = 16 

s _ 1 
l ~ 8 

s _ 3 
l 16 

s _ J _ 
l 32 

0 

0.1992 

0.1*131* 

0.71*69 

1.0397 

1/8 

.1708 

.3339 

.1*679 

.3969 

1/1* 

.1331* 

.21*62 

.271*8 

.1836 

1/2 

.0731 

.1198 

.1028 

.0397 

3A 

.0361 

.031*9 

.01*23 

.0236 


.0172 

.0233 

.0183 

.0101 


.0008 

.0012 

.0008 

.0003 
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TABLE 17.- NORMAL STRESSES oj FOR ALL 
ELLIPSES FROM (2X n b) Q - Concluded 


(d) a/b = 0.^000. 


X 

2b 

s _ 1 
1 16 

s _ 1 
l 8 

s _ JL - 
l 16 

s = JL 
7 32 

mm 

0.1007 

0.2031 

0.2700 

0,2214 

Sm S' 

• 0804 

.1503 

.1633 

.1088 


.0614 

.1072 

-.1013 

.0617 


.0331 

..052J 

. .0420 

.0238 

mfimm 

.0166 

..0248 

' .0186 

.0102 

wmmmm 

.0083 

'.0119 

.0086 

.0047 

■EH 

.0005 

.0007 

.0005 

.0003 



(e) a/b = O'. 8660. 


X' 

' s _ 1 

s _ 1 

! = J_ 

s = JL 

2b 

l 16 

1 8 

1 16 

z 32 

0 

0.0244 

0.0383 

■ 

0.0309 

0.0175 

1/8 

.0181; 

.0281 

.0218 

.0122 

lA 

.0137 

.0205 

.0155 

.0086 

1/2 

■ .0075 

.0109 

.0079 

.0043 

3 A 

' .0040 

.0058 

.0042 

.0023 

1.0 

.0022 

.0031 

.0022 

.0012 

2.0 

.0002 

.0002 

.0002 

.0001 
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TABLE 18.- SHEAR STRESSES t/t 0 FOR ALL 


ELLIPSES FROM (2X n b) 0 
(a) a/b = 0.0871 
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TABLE 18.- SHEAR STRESSES t/t 0 FOR ALL 
ELLIPSES FROM (2\ n b) 0 - Concluded 
(d) ' a/b = 0.5000. 




























NACA TN 2600 


TABLE 19.- COORDINATES OF NACA 63]_-012 AIRFOIL 
SECTION' TORSION BOX 




— = C 1 + 

7. 



cos 


2ruis 

1 


computed from 
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TABLE 20.- SERIES COEFFICIENTS FOR 
NACA 63i~ 01 2 AIRFOIL SECTION 


n 

V 

B n 

. ( 2 v>) 

1 

0.00288 

-0.00660 

2.31 

2 

.031*1*5 

-.011886 

1;. 86 

3 

-.00205 ' 

.00206 

7.38 

it 

.03068 

-.01870 

9.1;1 

5 

-.00386 - 

.00190 

11.79 

6 

.02553 

- -.01033 

, lli.10 

7 

-.001*09 

.00160 

16.90 

8 

.0205U , 

-.00598 

, 18.61 

9- 

-.00251* 

.00053 

21.31 

10 

.01613 

. -.00367 

23.12 

11 

-.0001*1* 

.00030 

26.1*6 

12 

.01221 

-.00237 

27.92 

13 

.00185 

— .0001*5 - 

32.25 

lit 

.00588. 

-.00131; 

33.79 

15 

.00377 

-.00061 

35.26 

16 

.00351 

-.00056 

38.06 

17 

.0051*9 

-.OOO93 

1*0.98 

18 

.00183 ' 

-.00006 

1*3.03 

19 

.00668 

-.00006 

1*1*. 76 

20 

-.00038 

.00015 

1*8.01 
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TABLE 21.- ANGLES OF TWIST q/{Y[/1 q G) AND 0/(M/lG) 
FOR NACA 63^012 AIRFOIL SECTION 


X 

0 

0 

2b 

MAoG 

M/IG 

0 

8.70 

1.00 

vu 

11.6k 

1.3U 

1/2 

12.33 

1*1*2 

1.0 

12.61 

i*U5 

1.5 

12.67 

ltlj.6 

2.0 

12.69 

l.lj.6 

U.o 

12.69 

1.U6 
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TABLE 22.- SUMMARY OF NORMAL STRESSES ct x ' FOR 
NACA 63i-012 AIRFOIL SECTION 


X 

2b 

S _ 1 
Z 16 

s _ 1 
1 8 

*~i|cn 

II 

s = _I 

Z 32 

s 1 

l h 

s _ 3 
l 8 

s = J_ 
Z 16 

1 = 15 

Z 32 

0 

- 1.90 

- 1.20 

- 0.83 

- 0.18 

-048 

0.83 

2.38 

5.13 

1A 

- 1.00 

-.82 

-47 

-40 

-.11 

.70 

1.07 

.72 

1/2 

-.36 

-.39 

-.25 

-.16 

-.07 

.32 

.35 

.22 

1.0 

-.09 

-.12 

-.09 

-.07 

-.03 

.08 

.07 

.Ol| 

1.5 

-.02 

-.03 

-.03 

-.02 

-.01 

.01 

.01 

.01 

2.0 

-.01 

-.01 

-.01 

-.01 

-.01 


0 

0 

l*.o 

0 

0 

0 

0 

0 


0 

0 


TABLE 23.- SUMMARY OF SHEAR STRESSES %/x Q FOR 
NACA 632-012 AIRFOIL SECTION 


X 

2 b 

B 

B 

s _ 1 
Z 8 

s _ 3 
Z 16 

s _ 7 
l 32 

s _ 1 
l h 

s _ 3 
Z 8 

np 

s _ 15 

Z 32 

s _ 1 
Z 2 

0 

9-85 

2.62 

1.90 

1.71 

1.63 

1.60 

1.71 

1.66 

1.90 

7.53 

lA 

748 

3.95 

2.26 

1.96 

2.23 

2.68 

3.77 

3.97 

. 01 * 

2.91 

1/2 

6.32 

1-.80 

3.18 

3.17 

3 .U 

3.21 

3.63 

1*.21 

.03 

3.10 

1.0 

6.19 

1 . 61 * 

3-30 

343 

3.38 

349 

3.82 

409 

-. 21 * 

2.92 

1.5 

6.4 

1.57 

3.31 

3.50 

346 

3.60 

3.85 

It . 05 

-.31 

2.89 

2.0 

6.13 

1.56 

3.30 

3.52 

348 

3.62 

3.87 

It . 05 

-.32 

2.89 

U.o 

6.12 

1.55 

3.29 

3.52 

349 

3.63 

3.88 

1*.06 

-.32 

2.8 9 












































(a) Values based on X^, and X^,. Comparison vitb values obtained by 
Karman and Chien (reference l) Is shown for £ - 0.10. 

Figure 6,- Variation of angles of tviBt 0(M/lG) vith x/2b for 
t angular boxes vith various 8 values. 
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Figure 8.- Variation of normal stresses -a x ' with s/l for rectangular 

boxes with. 0 « 0.05. 
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Figure 12.- Variation of normal atressi 
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Figure 19,- 

Variation 

of normal stresses 
"boxes with a/b «= 

* 

ct x ' with 
0.1738. 

s/l for elliptical 
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"boxes 
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0 ^' vfth b /1 for elliptical 
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Figure 25.- Variation of normal stresses a x ' with s/Z for 

MCA 63i-012 airfoil. 
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